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It is suggested that the study of multiplicity difference correlators between two well-separated 
bins in high-energy heavy-ion collisions can be used as a means to detect evidence of a quark-hadron 
phase transition. Analytical expressions for the scaled factorial moments of multiplicity difference 
distribution are obtained for small bin size with mean multiplicity s < 0.3 within Ginzburg-Landau 
description. It is shown that the scaling behaviors between the moments are still valid, though the 
behaviors of the moments with respect to the bin size are completely different from the so-called 
intermittency patterns. A universal exponent 7 is given to describe the dynamical fluctuations in 
the phase transition. 



One of the primary motivations of the study of high-energy heavy-ion collisions is to investigate the properties of 
quark-gluon system at extremely high temperature and high density. Such system may be in the state of quark-gluon- 
plasma (QGP), and with the expanding and cooling the system will undergo a quark-hadron phase transition and turn 
out to be hadrons detected in experiments. One of the theoretical aims is to find a signal about the phase transition. 
As is well-known for a long time, fluctuations are large for statistical systems near their critical points. Thus the 
study of fluctuations in the process might reveal some features for the phase transition. I 1 ' 2 1 Monte Carlo simulations 
I 3 1 on intermittency M without phase transition for pp collisions ^ show quantitatively different results on multiplicity 
fluctuations from theoretical predictions with the onset of phase transition.! 2 ] These different results stimulated a lot 
of theoretical works on multiplicity fluctuations with phase transition of second-order l 2,6 l and first-order I 7 ' 8 ! within 
Ginzburg-Landau model which is suitable for the study of phase transition for macroscopic systems. Most of these 
works give remarkable scaling behaviors between F q and F%, and there seems to exist a universal exponent v I 2 ' 6 - 8 !. It 
is suggested that the exponent v can be used as a useful diagnostic tool to detect the formation of QGP. In [7] In F q 
are expanded as power series of (5 1 / 3 , and it is shown that the coefficient of S 1 ^ 3 term can be used as a criterion for 
the onset and the order of the phase transition. All those works show the violation of intermittency patterns in the 
phase transition. 

It is known for a long time that the investigation of multiplicity fluctuations is very different in heavy-ion collisions, 
though the power-law dependence of hiF q on 8, F q oc 8~ Vq , has been found ubiquitous in hadronic an leptonic 
processes. ^ The main differences between heavy-ion physics and hadronic & leptonic ones on multiplicity fluctuations 
were noticed earlier in Ref. [10]. In Ref. [11] an alternative way was proposed to study the fluctuations by means 
of factorial moments of the multiplicity difference (FMMD) between two well-separated bins. This alternative is a 
hybrid of the usual factorial correlators M and wavelets I 12 ' because Wjk in Haar wavelet analysis is just the difference 
of multiplicities in two nearest bins. Present discussions, of course, will not be limited in the nearest bins. Let the 
two bins, each of size 8 2 and separated by A, have multiplicities ri\ and ri2, and define their multiplicity difference 
m = \n\ — tt-2 I ■ Scaled FMMD are defined as 

Fq = f q /fi, f q = 22m{m-l)---(m-q + l)Q m , (1) 

m 

with Q m the distribution of multiplicity difference which may be dependent on A, 8 and details of the process. 
Moments defined above are similar to but not the same as the Bialas-Peschanski correlators M F qiq2 , for T q may 
depend on both A and 6 while F qiq2 depends only on A. 

In Ref. [11], T q are numerically studied within Ginzburg-Landau model. The scaling behaviors between F q and 
F2, T q oc T2 , are shown with (3 q = (q— l) 7 and a universal exponent 7=1.099. 

In this paper, T q are studied analytically for very small bin size 8. Then the dynamical fluctuation components 

T^ dyn) of FMMD are defined. It is shown that both ]nT q and In J^ dyn) increase linearly with the bin size 8 when 8 
is very small, completely different from the usual intermittency behaviors of In F q which increase with the decrease 

of bin size. But the scaling laws between T q and T2, and between T q dv1 ^ and T^^ are still valid, although the 

corresponding j3 q and (3 q Ayri ^ are different. A universal exponent 7 for /3q dyn ^ is given which has no dependence on any 
parameter in the model and is different from that in [11]. 
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As an starting point, let us first discuss the trivial and simplest case. Suppose that the two bins considered are 
well-separated so that there is no correlations between them. Let the mean multiplicities in each bin are Si,S2, 
respectively. If there is no dynamical reason, the multiplicity distribution for each bin is a Poisson one 
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■exp(-Si) (i = l,2) 



From this distribution, one can deduce the multiplicity difference distribution as 

Try g -I 

P m {si, s 2 ) = cosh(- In -±)I m (2^^)e- Sl - S2 (2 - 

2 S2 

where I m (z) is the modified Bcssel function of order to, 
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FMMD for pure statistical fluctuations are 
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For simplicity, let us discuss the case with s\ = s 2 - This condition can always be satisfied if one chooses the two 
bins properly Since only to > q contribute to f q , the summation over to in last equation can be extended to m = 0. 
This summation converges very slowly because J m (2s) decreases with to approximately as s m /m\ for large to and 
small s, but the product to (to — 1) • • • (to — q+1) increases with to quickly. So contributions from all to > q must be 
taken into account. This will cause some difficulties in numerical calculations if one starts directly from the definition 
of the moments. In this paper, we will alternatively sum over to analytically, and then do numerical calculations 
from the final expression. In this approach, one can control the precision more easily in calculation. To complete the 
summation, one can introduce a generating function 



G(x,s) = 2e- 2s J2 x m Im(2s) , G q (x,a) = 



m=0 



diG(x,s) 
dxi 



With this function, ^ stat ) can be rewritten as 

/f at ) =G q (l,s) = G q (s) . 

Direct algebra shows that 

d 1 1 — exp(— xs) 



G(x,s) = 2e (x ~ 2)s 
with a° = {-l) l s 2i /{i\) 2 for i = 0, 1, • • and that 

j(stat) _ 2 e ~ s 



1=1 



dx l 
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where a\ can be calculated by recurrence relation from a?, a\ = sa^ , a\ = sd 



sa 



(8) 



r+d 1 .(*>2). 



Then one can get two specially important coefficients ak 



s q+2 rp^g mog ^. j m p 0r tant advantage of such 



calculations is that these formalisms facilitate analytical calculations for quantities in the range of very small bin size 
in which we are now interested. We will discuss it later in this paper. 

Now, we begin to discuss the FMMD in second-order quark-hadron phase transition. We use the Ginzburg-Landau 
description to specify the probability that s hadrons are created in the two dimensional, such as SrjSip, area 5 2 . In 
this description, the distribution of multiplicity is no longer a Poisson one and that for multiplicity difference is given 
by! 11 ! 
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Q m (S,r) = Z- 1 J V4>P m (5 2 r | <j> \ 2 )e- p M , (9) 

where r is an indication of lifetime of the whole parton system, T><p = nd \ <fi | 2 , Z = J T>(j>e~ F ^ and the free energy 
F \4>] = J S 2 dz [a \ <f> \ 2 +b \ <j) | 4 +c | d<f)/dz | 2 ] . 

As has been pointed out in [2, 6] that for small bin the gradient term in F[<j)] does not have any significant effect 
on the multiplicity fluctuation, so one can set c = 0. This setting means that 4> can be regarded as a constant over 
the area 5 2 . Of course, this is approximately true only when S 2 is very small. 

Substituting Q m (<5, r) into Eq. (1), one gets 



/,= J duG q (Txu)e xu - u J j due xu - u (10) 
with x ~\ a | 5/b related with the bin width S. Define' 6 ' 

poo 

J q {a) = duu q e au - u2 (11) 



which satisfies recurrence relation J q {a) = j.J q -i(a) + ^^-^-2(0:) and can be directly integrated for q — and 1, 
Ma) = ^e Q2 / 4 (l + erf(f )), Ji(a) = \ + f J (a). With J q (a), f q can be expressed as 

CO 

/, = J^i^zZbl^TM-ir-m (12) 

i=q 

with b\ constants, b q q = 1, especially. Notice that the second nonzero bf for fixed q is for i = q + 4. One can check this 
from the expression for G{x, s) and the recurrence relations for af. 

The scaled FMMD T q defined do contain contributions from statistical fluctuations, contrary to the usual scaled 
factorial ones. As a way to seek for the dynamical fluctuations, one can define the dynamical scaled FMMD as 

^ dyn) - Jfe • (13) 

J- q 

To make the definition sense, one should ensure that the mean multiplicity is the same for all the calculation of the 
moments concerned. In Ginzburg-Landau model, the mean multiplicity is s = tx Ji(x) / Jq(x). Then deviations of 

4 dyn) from one should indicate the existence of dynamical fluctuations. The three classes of moments defined in this 
paper can all be calculated directly. 

Up to now, all of moments are expressed as infinite sums and are exact within the model. The infinite summation 
will hinder us from an explicit formalism for interesting quantities. Now we focus on the range of very small bin 
size. As has been shown, the smallness of the bin size 5 is for the need of self-consistence of Ginzburg-Landau model 
adopted in this paper, otherwise the gradient term plays a role and cannot be set to zero. Experimentally, the bin 
size S can be chosen very small indeed. For example, experimental data' 13 ' show that the number of total produced 
charged particles is about 70 within a rapidity range about 7 in 200 A GeV S+Em collisions. The rapidity resolution 
in EMU01 experiments can be high up to 0.01. In two dimensional analysis as in this paper, the area S 2 considered 
can be so small that in that area the mean multiplicity satisfies s« 1. The mean multiplicity in single bin can still 
be much less than 1 even for Pb-Pb collisions in which the number of produced particles can reach 1500 or more. 
Because of such experimental facts, we can discuss only the cases with s <C 1 in the following, and our results can be 
checked directly in experiments. One can see that this condition will enable us to reach simple expressions for all the 
moments. 

For the pure statistical fluctuation case, terms excepts the leading term in Eq. (8) can be neglected, and one can 
easily get 

lnF^ = (q- l)(s-ln2) . (14) 

One can check that the relative contribution from all non-leading terms is about 1% for s = 0.3. 

For the moments with the onset of phase transition, it is a little complicated because of the integration in Eq. (10) 
over the whole range of s. But, one can see that the leading term in G q (s) plays a dominated role. One needs to 
notice that integrating u n term is associated with product of two factors (rx) n and exp(— t'xu — u 2 ), r' = r — 1. For 
small xt the first factor strongly suppresses the contribution. For larger xt the term exp(— t'xu) over-depresses the 
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contribution from the former. In fact, numerical results show that (sr) 4 J q+ 4(— t'x)/ J q {— t'x) is always of the order 
10~ 4 for xt < 0.5, corresponding to s ~ 0.3 for r = 10.0. So that the results will not be affected practically if only 
the leading term are kept for the calculation of the moments in small x region. Then to a good approximation, 

In T q = (q - 1) In + In J q (- T 'x) - q In Ji(-rV) , (15) 

\nT^ = (q-l)ln^^-+\nJ q (-T'x)-q\nJ 1 (-T'x) . (16) 

The behaviors of In T q and In T q Ayn ^ as functions of x from 0.005 to 0.05 are shown in Fig. 1 for r=2.0 and 10.0. 
The x range is chosen from the requirement s < 0.3 for t=10.0. One can see that both In T q and In J r g dyn - ) have linear 
dependence on bin size x. This dependence is completely different from the usual intermittency behaviors. This result 
can also be seen directly from last expressions for the moments if one substitutes J q (a) with \ (r(^-) + aT(^-)) 
for very small a. In small x approximation, 



In T q = const + 



=^17-9^1 x + Otf), (17) 



\nT<f yn) = const + t' 



r(l) 



r(i) r(2+i) 




.t + C>(:^) . (18) 



Numerical results show trivial scaling behaviors for In T q vs lnJF 2 and In ^Fq dyn ^ vs In T^ yn ^ in Fig. 2. Though hiJ- q 
have different ranges of values for different r, the scaling behaviors seem independent of the lifetime of the system. 
One can see weak dependence on r for (3 q from last equations. P q dyn ^ do not depend on any parameter in the model 
because the r dependencies in the local slopes are cancelled miraculously with each other in small x limit. More 
interestingly, f3 q dyn ^ can be well fitted by 

f3^ = (q-iy (19) 

with 7=1.3424, as shown in Fig. 3. But (3 q do not obey the same scaling law, as shown in Fig. 3 for the case with 
t = 10.0. The universal exponent 7 is different from that in [11]. But the difference does not mean any contradiction 
between present paper and [11], because they correspond to different quantities. The difference also comes from 
the different x regions discussed since (3 q depend on the fitting range. In this paper, the exponent 7 is completely 
determined by the general features but does not depend on any parameter of Ginzburg-Landau model used to describe 
the phase transition. The exponent 7 given here is very close to the exponent v given in former studies on multiplicity 
fluctuations for second-order phase transition. The slight difference between them comes from the different regions 
concerned. As shown in Fig. 3 of the first paper in [6], (3 q take minima at a = (x/2) 5 ~ 1 and increase with 
the decrease of a. Our result corresponds to a = 0. Since both v and 7 describe dynamical fluctuations in phase 
transition, they should be equal, as physically demanded. For experiments with s in single bin a little larger than 
0.3 the experimentally obtained 7 should be close to but less than 1.3424. Thus if experiments observe a scaling 
exponent 7 about 1.34 in a high resolution analysis, the onset of a second-order quark- hadron phase transition can 
be pronounced. 

It should be pointed out that even without dropping off non-leading terms, the exponent 7 in this paper will not be 
changed, because all those terms are related to higher orders of x and have no contribution to 7 which is connected 
with properties of the moments in the limit x — > 0. In this sense, the exponent 7 given here is exact and truly 
universal. 

In summary, scaled FMMD are studied analytically within Ginzburg-Landau model in a kinetical region with mean 
multiplicity in single bin less than 0.3 for second-order quark-hadron phase transition. The dynamical fluctuations in 
FMMD are extracted, which give the same physical contents as the usual scaled factorial moments. Scaling behaviors 
between scaled FMMD are shown, and a truly universal exponent is given. 
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Figure Captions 

Fig. 1 Dependences of In T q and XviTq on the bin width x for r = 10.0 and r = 2.0. 

Fig. 2 Scaling behaviors of In T q , In !Fq n ' vs In J 7 ^, In for the same choices of lifetime as in Fig.l. 
Fig. 3 Scaling exponent \n[3 q vs ln(g — 1). 
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